Dynamics of immersed molecules in superfluids 
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The dynamics of a molecule immersed in a superfluid medium are considered. Results are derived 
using a classical hydrodynamic approach followed by canonical quantization. The classical model, a 
rigid body immersed in incompressible fluid, permits a thorough analysis; its effective Hamiltonian 
generalizes the usual rigid-rotor Hamiltonian. In contrast to the free rigid rotor, the immersed body 
is shown to have chaotic dynamics. Quantization of the classical model leads to new and experi- 
mentally verifiable features. It is shown, for instance, that chiral molecules can behave as "quantum 
propellers" : the rotational-translational coupling induced by the superfluid leads to a nonzero linear 
momentum in the ground state. Hydrogen peroxide is a strong candidate for experimental detection 
of this effect. The signature is a characteristic splitting of rotational absorption lines. The loi — > lio 
line in hydrogen peroxide, for example, is predicted to split into three lines separated by as much 
as 0.01 cm -1 , which is about the experimental linewidth. 



I. INTRODUCTION 

The dynamics of a molecule moving within a bulk su- 
perfluid medium are not expected to be too different from 
its dynamics in vacuum. Like the vacuum, an unbounded 
fluid is homogeneous and invariant with respect to con- 
tinuous translations and rotations. Moreover, in the ex- 
otic case of superfluids, the very low level of excitation 
achieved at low temperatures implies that there is a van- 
ishing set of modes with which the immersed molecule 
can interact. These expectations are borne out in the 
spectroscopy of single impurity molecules in He nan- 
odroplets, where the superfluid medium has little effect 
on the rotational spectrum beyond modifying the molec- 
ular moments of inertia [2-4]. 

On the other hand, a low temperature superfluid 
medium is distinct from the nonrelativistic vacuum in 
that it lacks a basic symmetry: Galilean invariance. The 
superfluid condensate defines a preferred rest frame, an 
"ether" relative to which molecular velocities should be 
measured. This implies, in particular, that the energy- 
momentum relationship of molecules immersed in super- 
fluids can deviate from a parabolic form. Of more interest 
spectroscopically, the presence of a preferred rest frame 
also gives rise to the possibility of couplings between lin- 
ear and angular momenta that are strictly forbidden in 
an environment with Galilean invariance. Evidence of 
such couplings would provide one of the most direct ex- 
perimental signatures of the superfluid medium. 

This paper develops a simple and general theoretical 
framework for the rigid-body dynamics of a molecule im- 
mersed in a bosonic superfluid. The classical analysis of 
an immersed body, presented in Section II, shows that 
an incompressible fluid medium adds 11 parameters to 
the usual 4 parameters characterizing a free rigid body. 
As a result, an immersed body is typically chaotic. The 
classical model is quantized in Section III, and the quan- 
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turn Hamiltonian is discussed. The most interesting new 
feature is a term of the form p • J which directly couples 
the linear momentum of the molecule to its angular mo- 
mentum. This coupling requires chirality in the molecule, 
and corresponds intuitively to the action of a propeller. It 
fully lifts the usual (2 J+ l)-fold degeneracy of rigid-rotor 
states with angular momentum J, providing a clear spec- 
troscopic signature of the superfluid-induced rotational- 
translational coupling. An unusual consequence is that 
the ground state acquires nonzero linear momentum. In 
Section IV, the magnitude of the rotational-translational 
coupling is estimated for the chiral molecules HOOH and 
FOOF. The predicted splitting of rotational absorption 
lines, at least for HOOH, is comparable to the resolution 
of current experiments. Finally, in Section V, the results 
are discussed in the context of related work. 



II. CLASSICAL MECHANICS OF AN 
IMMERSED RIGID BODY 

We consider a rigid body immersed in unbounded, in- 
compressible fluid of uniform density pj. The body occu- 
pies a bounded, simply-connected volume V from which 
the fluid is excluded; the surface of this volume is denoted 
by dV. The velocity of the fluid is described by a poten- 
tial 4>(x), such that V/ = — V0, which satisfies V 2 = 
inside the fluid volume and |V0(x)| -> as |x| — * oo. 
On the moving surface dV, the fluid velocity also must 
satisfy (v/ — v s ) ■ n = 0, where n is the surface normal 
and v s is the surface velocity. Our goal is to obtain a 
classical Hamiltonian for the immersed body which can 
be readily quantized. The classical problem is treated in 
many hydrodynamics texts (e.g., Ref. [5]), to which we 
refer the reader desiring greater detail. 

We first choose an origin and coordinate axes fixed in 
the body; this defines the body frame. The origin of the 
body frame is at y in the lab frame, and the body axes are 
rotated by 7Z from the lab axes. The three components of 
y and the three Euler angles needed to parametrize the 
matrix 1Z completely specify the position of the body. 
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When y = and 1Z = 1, the body and lab frames coin- 
cide, and the body is said to be in the reference position. 
Let E t be the Euclidean transformation which moves the 
body from the reference position to its position at time 
t. Explicitly, 



E t : x ^ y(i) + n(t)x. 



(1) 



The time-dependence of V is simply Vt — E t Vo, where Vo 
is a fixed volume; similarly, dVt — E t dVo. The lab- frame 
velocity of a point E t x, which is fixed in the body, is 



-CEtx)=y + wx(ftx), 



(2) 



where u> is the body's angular velocity. The velocity 
potential at a point x which is fixed in the body frame, 
</>'(x) = (f)(E t x), satisfies fairly simple equations: V 2 ^' = 
outside Vo, |V0'(x)| — > as |x| — > oo, and 

V0'(x) -n = - (£ T v s (£ t x)) n 

'ft T y + (K T u>) x x) • n (3) 



for x on the fixed surface 8Vq, where the surface veloc- 
ity was calculated using Eq. (2). Indeed, the equations 
for <fi' are linear^ with inhomogeneous boundary condi- 
tions linear in lt T y and lZ T u>. The velocity potential 
can therefore be expressed as a linear combination, 

0'(x) = -(TZ T y)M^) - (£ T "to(x), (4) 
where the harmonic functions an( l satisfy 

VVv( x ) • n = e M • n, (5a) 
Vx M (x) • n = (e M x x) • n (5b) 

on dVo- These six functions characterize the fluid re- 
sponse to the six independent motions of the body; they 
depend only on the shape of dVo- 

For a free rigid body, the Lagrangian is just the kinetic 
energy. When the body is immersed, this is augmented 
by the kinetic energy of the fluid: 



6L 



P/|V0(x)| 5 



X/4 

J £x\l> t \V<l/{x)\*, (6) 



where the second integral is over the fluid volume V§, 
the complement of the body volume Vq. This is easy to 
evaluate using Eq. (4). The total Lagrangian, including 
the fluid contribution [Eq. (6)], then has the form 



L = \ (y T TZ ^TZ) 



M G 
G T I 



K T y 
K T u> 



as a sum of rigid-body and fluid tensors: M = Rd^+SM 
and so on. 

The rigid-body tensors are 

= J d 3 xp(x)5^, (8a) 

G$ = J d 3 xp(x)e^ x x Xl (8b) 

1$ = / d 3 xp(x) (|x| 2 ^ - Xli x v ) , (8c) 

where p(x) is the body-frame mass density of the body. 
The expressions for and J(°) are the usual ones, and 
G^ is typically made to vanish by choosing the body's 
center of mass as the origin of the body frame. The rigid- 
body Lagrangian then takes the familiar form 

L (0) = i y r Af (°) y + l -u T Ki^K T u. (9) 

The fluid tensors are 

SM^ = I d 3 xp f Vip fl (x)-Vip v (x), (10a) 

SG^ = f d 3 xp f V^(x) ■ V X „(x), (10b) 
SI„ V = / d 3 xp f Vx»{*)-VxM)- (10c) 

Note that the mass and inertia tensors M and I are sym- 
metric, while the rotational-translational coupling G has 
no obvious symmetry. In the most general case there 
are 6 + 6 + 9 = 21 parameters in the Lagrangian. We 
can place the body origin at a conveniently chosen point 
in the body, reducing this to 18. We can then rotate 
the body axes with respect to the body, leaving 15. By 
comparison, the free rigid rotor has 4 parameters. 

The classical Hamiltonian follows from the Lagrangian 
[Eq. (7)]: 



H = \ ( P ^ G t 



(11) 



(7) 



where p is the linear momentum, canonically conjugate 
to y, and J is the angular momentum in the lab frame. 
The components of J have Poisson brackets [J M , J„] cl = 
£nv\J\- It is convenient also to work with J' = TZ 7 J 
and p' = TZ T p, the angular and linear momenta in the 
body frame, in terms of which the Hamiltonian is par- 
ticularly simple. These variables have Poisson brackets 

[ J '^ J 'u] c i = -t^xJ'x and [ J ^KL = - e ^xP'x- Thc 
classical Poisson brackets relate to the quantum commu- 
tators in the usual way, as [A, B] = i[A, B] c \. 

To express Eq. (11) in a form more amenable to quan- 
tization, we first define new tensors a, (3, and 7 through 



The mass tensor M, the rotational-translational coupling 
tensor G, and the inertia tensor I can each be expressed 



a (3 
F 7 



M G 

g t i 



(12) 
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In the following section it is shown that the model can 
be quantized without complication provided that is 
symmetric. We can ensure that is in fact symmetric 
by choosing the body frame's origin correctly. For a free 
rigid body, the body's center of mass is the correct choice, 
but this is not true in general. The tensors a and 
transform in a particular way when the body origin is 
translated within the body (7 is unchanged); we can use 
the transformation law to find the appropriate choice of 
origin. Specifically, consider shifting the body origin by 
d, so that all body- frame coordinates transform asxH 
x — d. Then a direct calculation from Eqs. (8) shows that 
the rigid-body tensors transform as 



/M(0) £(o)\ 



i 0\ /M(°) 6<°>\ (t -ed 
ed i) [d^ T /(°) )\0 1 



where the antisymmetric matrix ed has matrix elements 
(ed) M „ = e^ u \d\. It is not hard to see that the fluid 
tensors transform in exactly the same way. Inverting 
the relation, we find that 0^0 + (ed)7- Now, the 
vector b with components b\ — e^xP^u is the zero vector 
exactly when is symmetric, and it transforms asb^ 
b+ { (Try)! — 7} d. The matrix multiplying d is positive 
definite, because 7 is, so it can be inverted. The result is 
that if is not symmetric, it can be made so by shifting 
the origin of the body frame by d = — {(Try)! — 7} b 
within the body. We will assume that this has been done. 

With the body origin fixed, rotation of the body axes 
can be used to diagonalizc 7, leaving it in the form 
7 = diag(2A, 2B, 2(7). (Since transforms as a ten- 
sor under rotation, it will remain symmetric.) Note that 
A, B, and C are effective rotational constants, shifted 
from their rigid-body values by the body's interaction 
with the fluid. Finally, the symmetric matrices a and 
can be broken down into their scalar and rank-2 spherical 
tensor components: 



a = ai 0) i + 



(14) 



q=-2 



and analogously for 0, where 
1 



-1 



M { 2\ = ^ I ±i -1 1. 



1 ±i 0> 




l ±2 



(15) 



0, 



These matrices are defined to have nice rotational prop- 
erties; in particular, 



(16) 



(2) 

where the T> pq are rotation matrices, given by known 
functions of the Euler angles [6]. The Hamiltonian then 
has the form 



H 



1 



«oV + 0o 'p ■ J + n rr {A, B, C) 



?(o) 



^ 2 

q 



U^ + J20 q 2) e q 2) , (i7) 



where 7i rr is the rigid rotor Hamiltonian, 

H rr (A, B, C) = AJ' X 2 + BJ l2 + CJ' Z \ 



(18) 



and n( 2 ) and 0^ are spherical tensors of rank 2, with 
components 

nf) = (p'f^V; ew = ( P ') T ^fJ'. (19) 

More precisely, II^ 2 ) and O^ 2 ) transform as spherical ten- 
sors of rank 2 under rotations generated by J'. 

Before proceeding to the quantum case, we consider 
the classical dynamics of the model. The Lagrangian is 
invariant under arbitrary translations (y y + a) and 
system rotations (y 1— > Oy, TZ OTZ), but not under 
body rotations (7Z OTZ) alone. Therefore p and the 
total (spin plus orbital) angular momentum J + y x p 
are constants of the motion, but J is not, in contrast 
to the free rigid rotor. We can use these symmetries 
to simplify the system. Four constants of the motion, 
which have vanishing Poisson brackets with Tt and with 
each other, are p • J and the components of p. By fixing 
these constants and eliminating the coordinates which 
are conjugate to them (the components of y and one 
of the three Euler angles), we restrict the Hamiltonian 
to a reduced phase space; the equations of motion for 
the remaining degrees of freedom are unchanged by the 
reduction. In this case, the original six degrees of free- 
dom are reduced to two, so the reduced phase space is 
four-dimensional. It can be parametrized by the direc- 
tions of p' and J', expressed in polar coordinates by the 
four angles (9 p ,<f) p ,8j,(f)j). The energy shells of the re- 
duced phase space (i.e., the level sets of Ti.) are three- 
dimensional. 

The integrability of a Hamiltonian system with two 
degrees of freedom can be ascertained using a Poincare 
section, essentially by inspection. To illustrate that the 
immersed rigid body is generally chaotic, we display a 
Poincare section of the reduced phase space in Fig. 1. 
The figure was generated by repeatedly integrating the 
equations of motion, starting at many different points 
within the TL = 2 energy shell, and plotting the point 
(<j)j,cos6j) each time a trajectory crossed the plane 
9 p = it/2. The Hamiltonian parameters were a = 1, 
— diag(l, 1, — 2), and 7 = diag(l, 2, V5); the con- 
served quantities were set to p • J = 0.1 and p = 0.1e z . 
If the system were integrable, each phase space trajec- 
tory would be confined to a two-dimensional torus, and 
would appear in the figure as a densely dotted curve. 
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FIG. 1: Poincare section of the reduced phase space of an immersed body, showing both chaotic and quasi-integrable regions. 



Instead, the figure shows that the phase space contains 
chaotic regions, where trajectories ergodically visit three- 
dimensional volumes, as well as quasi-integrable regions, 
where trajectories are confined to two-dimensional tori. 
It should be noted that this set of parameter values cor- 
responds to very strong coupling between linear and an- 
gular momenta, and was chosen to highlight the presence 



of chaos. The physical cases we will consider in this pa- 
per are much more nearly integrable. Nevertheless, the 
existence of classical chaos suggests that the quantum 
spectrum, unlike that of the free rigid rotor, cannot be 
determined algebraically. We will need to use an approx- 
imate method, in this case a perturbative expansion in 
powers of the linear momentum. 



III. QUANTUM HAMILTONIAN 

The quantization of the classical Hamiltonian 
[Eq. (17)] is straightforward, and proceeds as in any 
standard angular momentum text (sec, for instance, 
Ref. [6]). There is only one subtlety. Since p' and J' 
have a nonvanishing Poisson bracket, their quantum 
counterparts do not commute, and the quantization 
of terms like p'^J^ suffers from an ordering ambiguity. 
We can sidestep this problem because we only need 
to quantize the components of 0( 2 ); since [^f^j^ is 

" (2) 

^-antisymmetric, while each matrix M q is symmetric, 

(2) 

the operators Q q are independent of the relative 
ordering of p' and J'. For the remainder of the paper, p, 
H, etc., will refer to the appropriate quantum operators 
(or their eigenvalues) unless otherwise noted. 

Breaking the Hamiltonian down into H.q and 57i, re- 
spectively the first three and last two terms of Eq. (17), 
allows us to proceed perturbatively. Since 5H vanishes 



for a free rigid rotor, this separates out the effects of the 
fluid. The p • J term also is absent for a free rotor, but 
this term commutes with p 2 and H rr , so it can be treated 
exactly. 

Because p is conserved, we can fix p = ke z , k > 0, with 
no loss of generality. A convenient basis for the remain- 
ing (angular) degrees of freedom is the set of simultane- 
ous eigenvectors of J 2 , J' z , and J z ; the eigenvector with 
respective eigenvalues J(J+ 1), K, and M is denoted by 
| JKM). The tensor operators [Eq. (19)] become 

Il q V = ^vi 2 q \ (20a) 

e?> = * + - \j + t> { S) < ( 20b ) 

where J± = J x ±iJ y . Their matrix elements are express- 
ible in terms of Wigner 3-j symbols [6] : 



■5 



(J 2 K 2 M 2 \U^ \J X K X M X ) = _fc 2 (-l)^+ M V(2J 2 + Wh + 1) 



J* 


2 


Jl ) 


J 2 


2 




\K 2 


Q 




\M 2 








(-1) J 



Kl ( J 2 2 Ji 
^ 2 g -Jfi 



J 2 M 2 



n (2) 



and 



J 2 2 Ji 
K 2 q -K x 



(J 2 K 2 M 2 \ 9< 2 > \J 1 K 1 M 1 ) = fc(-l)^ +M V(2J 2 + l)(2Ji + l) ( 

{lV(J2-M 2 )(J 2+ M 2 + l)( M J 2 l 1 \ \ 



-\^{J 2 -M 2 + l){J 2 + M 2 ) 



J 2 2 Ji 
M 2 - 1 1 -Mi 



= (-1) 



A'l 



J 2 2 Ji 
if 2 g -Xi 



J 2 M 2 



6 (2) 



JiMi 



(21a) 



2M 2 / J 2 2 Ji 



76 \M 2 -Mi 



(21b) 



The simple dependence of the matrix elements on Ki , K 2 
and q is a consequence of the tensorial nature of the oper- 
ators, and is guaranteed by the Wigncr-Eckart theorem. 
The reduced matrix elements vanish unless M 2 = Mi , by 
inspection, so J z commutes with both Ho and SH. This 
is to be expected from the classical conservation of p • J, 
which equals kJ z for our choice of p. 

The eigenvectors of Ho are the simultaneous eigenvec- 
tors of the rigid rotor Hamiltonian and J z . Since H rr 
can be written in a normalized form, 



Tt rr (A, B, C) 



l -(A + C)J 2 + l -{A-C)H r 



,(1,K,-1), 

(22) 

the eigenvectors depend only on the asymmetry parame- 
ter k = (-A + 2B - C)/(A - C). The eigenstates evolve 
continuously with k; they become eigenstates of J' x and 
J' z in the prolate (k — > —1) and oblate (k — ► 1) limits re- 
spectively, with corresponding eigenvalues K_i and K + i. 
Following standard notation, we denote these states as 
|rM), where r € {Jk^k +1 } = {0 00 , loi, In, lio, ■ • ■} 
indicates both the angular momentum and the two lim- 
iting values of K, and M is the eigenvalue of J z . The 
unperturbed energies are 



E (0) _ 



1 



+-(A-C)£ t (k) 



(23) 



written in terms of the normalized rigid rotor eigenvalues 
£ t (k). 

We are now in a position to find the perturbed energies 
to any order in SH. The expansion can be reordered as an 
expansion in the momentum k. The first-order correction 
comes from the diagonal elements of &^ 2 \ which turn out 
to be proportional to M for each r. A given level has 



energy 



E TM (k) = ^(A + C)J(J + 1) + ^(A-C)£ t (k) 

+B T Mk + X -A TM k 2 + 0(fc 3 ), (24) 



where B T and A t m are calculated from matrix elements 
of SH in the usual way. In particular, 



A T M 



4 0) + ^E4 2) K 2) l 



ttM 



(25a) 



k 2 



EWE 



©„7 



q' Jtt'M\^1 



rM 



q,q' 



,(o) 



(25b) 



where (Q) T t'm is shorthand for (tM \Q\t'M). Higher- 
order terms are suppressed by additional factors of kf3q 2 ^ 
or k 2 a^: for the cases we consider in the next section, 
these factors are less than 10~ 3 , so it is reasonable to 
truncate the expansion here. Clearly A~^ is an effective 
mass; it reduces to the mass of the body in the absence of 
the fluid. The constant B T is a pseudoscalar associated 
with the family of levels |rM). It is nonzero only in the 
presence of the fluid, and only when the immersed body 
is chiral. Judging from the spectrum, B T measures the 
tendency of the immersed body to have its linear and 
angular momenta aligned. When this tendency is strong, 
the body behaves as a kind of "quantum propeller." 



IV. MOLECULAR SURFACE MODEL AND 
NUMERICAL RESULTS 

In order to apply these results to a real system, we 
must specify the fluid density p/, model the surface 
dV and mass density p(x) which characterize the im- 
mersed body, and calculate the tensors M, G, and / 
from these inputs. In this section we consider the in- 
teresting case of an immersed molecule in superfluid he- 
lium. The fluid density pf is approximately 0.1 amu-A~ 3 . 
The mass density of the molecule can be represented as 
p(x) = X)i m i^( x — x i)i where the atomic coordinates 
Xi are known. Helium is kept away from the molecule 
primarily by the short-ranged Fermi repulsion between 
molecular electrons and helium electrons. This results 
in a smoothly varying helium density, interpolating from 
zero near the molecule to pf far away. Our model approx- 
imates this smooth variation by a discontinuous jump, 
localized on the imaginary surface 8Vq. The optimal 
choice of surface is not obvious, but some criteria are 
clearly important. The surface should share the sym- 
metries of the molecule; it should be smooth, since the 
electronic densities are smooth; and its size and shape 
should be physically reasonable: when compared to the 
actual helium density profile, the surface should approx- 
imate the surface p — \pf- With these criteria in mind, 
we model the surface dVo for an immersed molecule in 
the following way. First, we place a sphere of radius 
Ri at each atomic coordinate x^. Each radius is pro- 
portional to the van der Waals distance between a he- 
lium atom and an atom of type i: Ri — aRne-i, where 
the parameter c will be of order 1. (For our calcula- 
tions, we used i?Ho-H = 2.60 A, -Rho-f = 2.87 A, and 
Rue-O = 2.92 A [7].) The union of these spheres forms 
a cuspy volume. The surface radius of the body volume, 
relative to an origin within the molecule, is represented 
in polar coordinates as r s (9,(f)). We then smooth the 
function r s using the rotationally invariant linear oper- 
ator which maps Y/ m (6>, (f>) \— > exp (— al) Yi m (6, <f>), with 
a = 0.1, to yield the final surface dVo. (Many other 
smoothing operators would do as well; our results are 
not very sensitive to the choice.) This prescription gives 
a single-parameter family of smooth surfaces which share 
the symmetry of the molecule; the parameter c can be ad- 
justed to exclude superfluid from an appropriately sized 
volume. 

Once the surface is fixed, the fluid tensors can be cal- 
culated either numerically or analytically. We used a 
perturbative analytical method, valid for nearly spheri- 
cal surfaces, where the surface radius is a weakly varying 
function of polar angle: r s (6, 4>) = ro (1 + e(6, (f>)), with 
eCl. The fluid response functions an d Xm ( see Sec- 
tion II) can then be calculated as power series in e, as can 
the fluid contributions to the tensors. This calculation 
has been carried through to second order, which is the 
lowest order at which nontrivial rotational-translational 
coupling (i.e. a nonzero value for B T ) is seen. The cal- 
culation and results are described in the Appendix. 
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(0.10, 0.49, 1.53) 




FIG. 2: Structures of the chiral molecules HOOH and FOOF. 
Atomic coordinates (x, y, z) are shown, in A, for two atoms 
in each molecule. The undisplayed coordinates are obtained 
from these by a symmetry transformation, R v {k) : (x, y, z) 
{-x,y, -z). 



We present data for the low-lying levels of two different 
molecules: hydrogen peroxide (HOOH) and dioxygen di- 
fluoride (FOOF). These molecules are depicted in Fig. 2, 
with structural parameters taken from Rcf. [8] . Each has 
a Ci point group, which we have taken to be Ci(\j) by 
an appropriate rotation of body axes. The allowed ro- 
tational levels all have the same parity under rotation 
by 7T around the y-axis: either even (0oo, In, 2o2, ■ • •) or 
odd (loij liO; 2i2, • • •), depending on the symmetry of the 
joint electronic/nuclear wavefunction. In the case where 
allowed levels are odd, the J = state is forbidden, so the 
ground state of the immersed molecule has J — 1, and, 
by Eq. (24), a nonzero linear momentum. The ground 
state momentum is k = \B T0 / A TOi+1 \, where r is the 
ground state with J = 1. 

The results for HOOH and FOOF are presented in Ta- 
bles I and II respectively. (Energies are given in cm' 1 or 
K; these conversions are made using ks = he = 1, con- 
sistent with usual spectroscopic notation. Momenta are 
given in A -1 ; this conversion is made using K = 1.) The 
surfaces considered can be parametrized by their mean 
radius r , shown in the first column. The values dis- 
played correspond to the parameter range 0.5 < c < 1.5 
for each molecule. Although the inverse effective masses 
A t m are level-dependent, this dependence is very weak, 
and A t m ~ a °^ for the tabulated cases. This approx- 
imate value is displayed, as A, in the second column. 
The remaining columns list the zero-momentum energies 
E T (0) and chiral splitting constants B T for the C2(y)-odd 
levels t = loi and r = lio- 

To put the results in context, consider the physical 
parameters of current nanodroplet experiments, together 
with the hydrogen peroxide results. Ignoring finite-size 
effects, the thermal momentum distribution is expected 
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TABLE I: Results for HOOH. 



ro 


A 


^ioi(O) 




Ei w (0) 




1.86 


0.94 


1.64 


-2x10 


-4 


9.68 


3.2xl0~ a 


2.47 


0.89 


1.61 


-2x10 


-4 


8.94 


4.0xl0" 3 


3.06 


0.83 


1.57 


-3x10 


-4 


8.28 


6.2xl0" 3 


3.63 


0.76 


1.52 


-5x10 


-4 


7.73 


8.1 xlO" 3 


4.21 


0.67 


1.46 


-6x10 


-4 


7.35 


9.1xl0~ 3 


4.78 


0.59 


1.38 


-5x10 


-4 


7.13 


8.6xl0~ 3 


(A) (■ 


cm 1 -A' ! ) (cm x ) 


(cm" 1 - 


A) ( 


cm" 1 ) 


(cm" 1 -!) 



TABLE II: Results for FOOF. 
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to be Maxwcllian, i.e. dP/dk = Ck 2 cxp [-\Ak 2 /Tj . 
The nanodroplet temperature is T = 0.4 K. This corre- 
sponds to a typical thermal momentum of k rms = 1 A -1 
for HOOH. The ground-state momentum fco, on the 
other hand, is only about 0.01 A -1 , and is completely 
overwhelmed by thermal fluctuations. The r = Iqi and 



t = 1 io energy levels will each acquire fine structure by 
splitting into three evenly spaced sublevels with different 
values of M; we estimate the spacing for these two sets 
of sublevels to be 2-5 x 10~ 4 cm" 1 and 3-8 x 10" 3 cm" 1 
respectively. Moreover, the center level (M = 0) in each 
triplet is sharp, while the wings (M = ±1) are broadened 
by the thermal spread in momentum. Because the exper- 
imental linewidth for rovibrational transitions is small, 
on the order of 0.01 cm -1 , we expect this fine structure 
to be nearly resolvable in the loi — > lio absorption line 
for HOOH. If resolved, the line would appear as a sharp 
central peak, with broader peaks symmetrically placed 
on either side, at ±0.003 - 0.008 cm" 1 . 

More generally, the resolvable fine structure of the ab- 
sorption line for a dipole transition r a — ► r b will de- 
pend on the relative values of kB a , kB bl and the natural 
linewidth T. There may be three peaks (if kB b >T> 
kB b - kB a ), 2J a + 1 peaks (if kB b - kB a > T > kB b ), 
3(2J a + 1) peaks (if kB b and kB b - kB a > T), or more 
complicated possibilities when some peaks overlap. For 
the loi — > lio transition in hydrogen peroxide, our the- 
ory predicts that the initial level spacing is very small, 
while the final level spacing is approximately resolvable, 
for typical thermal values of k. The expected line shape 
for this scenario is schematically depicted in Figure 3. 
There are three peaks, corresponding to the three possi- 
ble final values of M. The central peak has the natural 
linewidth, since the M — energy has no first-order de- 
pendence on k, while the M = ±1 wings are thermally 
broadened. 



V. DISCUSSION 

In this paper we have presented a classical model for 
the behavior of a molecule immersed in superfluid; and 
we have shown that the quantized version of the model 
has interesting features which should be, under certain 
conditions, spectroscopically detectable. In particular, 
chiral molecules can act as "quantum propellers," cou- 
pling linear and angular momenta via their interaction 
with the superfluid medium, and this leads to a charac- 
teristic splitting of spectral lines. Our model simplifies 
the superfluid dynamics greatly, which presumably af- 
fects the quantitative accuracy of our results. We will 
conclude by discussing possible remedies for the short- 
comings of our model. 

One important physical effect that we have ignored 
is the formation of structured shells of helium atoms 
around a solvated impurity, which leads to a nonuni- 
form fluid density. The model can be naturally extended 
to include this, by allowing the fluid density pf to vary 
spatially, while maintaining a time-independent profile 
in the body frame. This approach has been used by 
other authors to estimate the superfluid-induced increase 



in the moments of inertia of linear or highly symmetric 
molecules [9, 11, 12]. The effect of a nonuniform density 
can be dramatic. For instance, SF 6 has I = 180 amu-A 2 
in vacuum, which increases by AI = 310 amu-A 2 when 
the molecule is immersed in a helium droplet. The model 
presented here gives A J < 25 amu-A 2 , an order of mag- 
nitude too small, confirming the assessment by other au- 
thors that a rotating ellipsoid model can seriously under- 
estimate AI [11, 12]. By contrast, a calculation using 
a nonuniform density predicts AI =170 amu-A 2 , which 
is more than half of the experimental value [11]. The 
comparison demonstrates that a rigid body with solva- 
tion shells attached drags more mass when it moves than 
one without them. It is tempting to speculate that the 
calculated chiral splitting constants B T will also increase 
when a nonuniform density is allowed; but this may not 
be the case, for two reasons. First, since the off-diagonal 
Hamiltonian matrix $ is equal to — M~ X GI~ X to lowest 
order in G, increases in the rotational-translational cou- 
pling tensor are counteracted by corresponding increases 
in the hydrodynamic mass and inertia tensors. Second, 
attaching solvation shells to the immersed body simply 
builds a larger rigid body, to some extent, and the con- 
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FIG. 3: Schematic representation of the predicted loi — ► lio absorption line in hydrogen peroxide. The strength of absorption 
is plotted as a function of the photon wavenumber, v. 



stants B T always vanish for a rigid body. Therefore, while 
allowing a nonuniform fluid density may alter our results, 
it is difficult to predict the outcome without performing 
a full calculation. 

By studying the dynamics of immersed molecules in 
bulk superfluid, we have ignored the finite-size effects as- 
sociated with immersion in a cluster or droplet. Our 
approach can be extended to handle reasonably large 
droplets with no fundamental difficulties. Lchmann has 
carried out a related calculation for the case of a linear 
molecule in a spherical droplet, taking into account the 
hydrodynamic effects which are present for an elliptical 
molecular surface (the modification of the rotational con- 
stants and the a matrix, in our notation), together with 
the effective interactions between the molecule and the 
droplet surface [10]. ft would be interesting to extend the 
current work along similar lines. We expect that such an 
extension would lead to results qualitatively similar to 
those presented here. Because parity is a good symme- 
try for a superfluid droplet, as it is for bulk superfluid, 
parity would still be broken by the addition of a chiral 
impurity molecule, and this would manifest itself in a 
characteristic splitting of the rotational spectrum. 

Microscopic details of the superfluid structure are 
also neglected by a continuum hydrodynamic approach. 
These details are better treated by path-integral and dif- 
fusion Monte Carlo methods, which have proven success- 
ful in predicting the rotational constants of immersed 
molecules [9]. It is likely that the hydrodynamic mass 
and rotational-translational tensors described here can 
also be extracted from such calculations, and this infor- 
mation would be a useful complement to the hydrody- 



namic results. Monte Carlo methods could also be used 
to address another neglected phenomenon: the quantum- 
mechanical tunnelling between left-handed and right- 
handed forms of a chiral molecule. This may play a 
pronounced role in hydrogen peroxide, where the tor- 
sional ground state (a symmetric superposition of left- 
and right-handed forms) is 11 cm -1 below the first (anti- 
symmetric) excited state [13]. Because the splitting is so 
large, comparable even to the rotational level spacing, the 
enantiomers are strongly mixed, making questionable our 
rigid-body treatment of the molecule. We expect the sur- 
rounding superfluid to suppress tunnelling, but a detailed 
calculation, without the assumption of rigidity, is needed 
for a quantitative assessment. (Some experimental sup- 
port is provided by the work of Nauta and Miller, who 
studied immersed HF dimers; they found that the bar- 
rier to interchange tunnelling was raised significantly by 
the presence of the superfluid [14].) The main difficulty 
in applying Monte Carlo methods would be finite-size ef- 
fects: since translational symmetry is strongly broken in 
small droplets, while rotational symmetry is preserved, 
the finite-size corrections to M and G are presumably 
more drastic than those to I. A number of other micro- 
scopic techniques have been applied to molecular impu- 
rities in superfluid droplets (see Ref. [9] for a recent re- 
view), such as density functional theory, and these more 
elaborate methods may also be applicable to the phe- 
nomena we have described. 

Finally, experimental data would be extraordinarily 
useful in refining the current model. The spectroscopic 
signatures we have described should be present for any 
chiral molecule, but not all candidate molecules will have 
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the large splitting constants necessary to resolve the 
fine structure. Nanodroplct experiments have been con- 
ducted using many different species of impurity molecule; 
Ref. [4] contains an exhaustive list. However, almost none 
of the studied molecules are chiral, so it is not surprising 
that no quantum propeller has yet been seen. We have 
suggested hydrogen peroxide (HOOH) as one strong can- 
didate. Ideally, the present paper will provide sufficient 
impetus for more detailed investigations, both theoretical 
and experimental. 



and so forth, where now all derivatives are evaluated at 
r = ro. Moreover, V 2 xp^ — and |V^^(x)| — > as 
|x| — ► oo for each n. The zeroth-order equation can be 

solved for ip$\ and in general the equation of order e n 
can be solved for ip^ once all the V^™^ with m < n are 
known. Therefore ip^ can be calculated to any order in 
e, though the process rapidly becomes tedious. Similarly, 
the rotational response function Xn satisfies the simpler 
equation 
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APPENDIX A: CALCULATION OF FLUID 
TENSORS 

For approximately spherical surfaces, it is possi- 
ble to calculate the hydrodynamic mass, rotational- 
translational coupling, and inertia tensors as perturba- 
tion scries in the deviation from sphericity. Consider a 
surface at r s (6, (p) = r (1 + e(6, <p)). The surface normal 
times the area element is 



nda 

rl<PVL 



(1 + e) (e r + i{e r x J)) (1 + e), 



where 



J = ie g 



1 d 



d 



sin 6 d</> ie<p d9 



(Al) 



(A2) 



is the usual vector of first-order differential operators on 
the sphere. For any scalar function $, 



nda 

r%d 2 n 



(l + eJ^ + ijS.Je, 



dr 



(A3) 



with derivatives evaluated at r = r s . Applying this to 
Eq. (5a) for the translational response function ip^ yields 

dtp 1 1 

~dV = n »- ^TT^ J ^ ' Je + TT7 ze - A ^ (Je)A ' 

(A4) 

where = e p • e r , and with derivatives evaluated at 
r = r s . If ipp is expanded in orders of e, as W = iffl + 
ipjP + . . ., and functions evaluated at r s are expanded 
in Taylor series around ro, then Eq. (A4) gives a single 
equation at each order in e: 



f i = --o(Je),-^ F Jx,-Je (A6) 



at r = r s , which follows from Eq. (5b); and this can be 
used to find Xn to anv order in e, in exactly the same 
way. 

Once the response functions are in hand, Eqs. (10) can 
be evaluated to obtain the fluid tensors. This, too, is 
done order- by-order in e. The results are best expressed 
as angular averages. Through second order in e, 



^Pf r o 




(A7a) 



= l^M-^^y (A7b) 



dr / ' 



(A7c) 



where (/) = j- J /(Q)c/ 2 f2, and functions are evaluated 

at r . Here the first-order fields ip^ and xft are deter- 
mined by their radial derivatives at ro, which are 



dr 
dr 



-i(Je) M . 



(A8) 



dr 



r e 



d^ 




dr 


= n M , 


9 2 4 0) 


1 


dr 2 





(A5a) Note that SG is antisymmetric to first order. Because 
the antisymmetric part of G is eliminated by the cor- 
j e v rect choice of body origin, it was necessary to carry the 

calculation of 5G to second order to obtain nontrivial 

(A5b) rotational-translational coupling. 
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